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We study the dynamics of two interacting Bose-Einstein 
condensates, by numerically solving two coupled Gross- 
Pitaevskii equations at zero temperature. We consider the 
case of a sudden transfer of atoms between two trapped states 
with different magnetic moments: the two condensates are ini- 
tially created with the same density profile, but are trapped 
into different magnetic potentials, whose minima are verti- 
cally displaced by a distance much larger than the initial 
size of both condensates. Then the two condensates begin 
to perform collective oscillations, undergoing a complex evolu- 
tion, characterized by collisions between the two condensates. 
We investigate the effects of their mutual interaction on the 
center-of-mass oscillations and on the time evolution of the 
aspect ratios. Our theoretical analysis provides a useful in- 
sight i nto the recent expe rimental observations by Maddaloni 
et al, |cond- mat/0003402. 
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I. INTRODUCTION 

Several papers have been recently devoted to the study 
of multiple-species interacting Bose-Einstein condensates 
(BEC) jl]-^. The first experiments involving the in- 
teraction between multiple condensates have been per- 
formed at JILA, with two condensates in the hypcrfinc 
levels \F = 2, m f = 2) and |1,-1) of 87 Rb confined 
in a Ioffe-type trap 0- Other experiments have been 
performed with the states |2, 1) and |1,— 1) in a TOP 
trap, in a situation of nearly complete spatial overlap of 
the two condensates. In fact they have almost the same 
magnetic moment and therefore experience two (slightly 
displaced) trapping potentials with the same frequencies 
The resulting dynamics reveals a complex struc- 
ture and it is characterized by a strongdamping of the 
relative motion of the two condensates H . These results 
have been analyzed in ||, by identifying two dynami- 
cal regimes: (i) a periodic motion with a slow frequency 
(with respect to the trap frequency), for small displace- 
ment between the traps; (ii) a strong non-linear mixing 
which leads to a damping of the relative motion, for larger 
displacements (but still smaller than the size of the ini- 
tial condensate). By comparing the numerical solution 
of the Gross-Pitaevskii equation with the JILA experi- 



ment ||, the authors of Ref. || find a fair agreement, 
although small oscillations remain undamped for longer 
times compared to the experimental observation (^j^]. 

A mixture of 87 Rb condensates in different mf states in 
a TOP trap has been experimentally investigated also by 
another group Q , but no effects of the mutual interaction 
have been observed. 

Recently we have reported some experiments per- 
formed at LENS with two interacting condensates [[[), 
which reveal several new features with respect to those 
analyzed in Refs. jjfl^]. In fact, we have considered a 
sudden transfer of atoms between the states |2,2) and 
|2, 1) of 87 Rb, having different magnetic moments. The 
1 2, 1), which is created with the same density distribu- 
tion and at the same position of the |2,2) condensate, 
feels a magnetic potential whose minimum is displaced 
along the vertical y axis by a distance much larger that 
the initial size of the condensate. Therefore, its dynam- 
ics is characterized by large center-of-mass oscillations. 
The two condensates periodically collide and these col- 
lisions have significant observable effects. Among them: 
i) the shape oscillations of |2, 2) condensate, triggered by 
the sudden transfer to the |2, 1) state, are significantly 
enhanced; ii) the center-of-mass oscillation frequency of 
the |2, 1) condensate is shifted upwards. 

By means of preliminary numerical simulations, we 
have already shown that the basic features of these exper- 
imental observations can be explained within the Gross- 
Pitaevskii (GP) theory (!]]. In this paper we better anal- 
yse the theoretical scheme and describe in more detail 
the numerical procedures and the approximations made. 

The paper is organized as follows: in Section O, we 
introduce the two coupled GP equations for the dynam- 
ics of the condensates and disc uss the approximations 
used to solve them. In Section III we show the results 



for the center-of-mass motion, stressing the effect of mu- 
tual interactions in the oscillation frequency. In Section 
|lV| we discuss the collective modes of each condensate 
and the time evolution of the aspect ratios. In Section 
we simulate the expansion of the two condensates af- 
ter switching-off the trapping potential. In this way, we 
can compare our theoretical results with the experimen- 
tal data |1| for the center-of-mass motion and the aspect 
ratios. 
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II. GROSS-PITAEVSKII THEORY FOR TWO 
COUPLED CONDENSATES 

At very low temperatures, when there is no thermal 
cloud, the evolution of the two interacting Bose-Einstein 
condensates can be described by two coupled Gross- 
Pitaevskii (GP) equations (i — 1,2): 



h z v 2 

2m 



V* 



E 

.7=1,2 



Airfi 2 a 



*<• (1) 



Here, \&i and ^ are the order parameters of the |2, 1) 
and |2, 2) condensates of 87 Rb (hereafter called |1) and 
1 2)). By solving Eq. (Q) we want to simulate the situ- 
ation in which a certain amount of atoms are suddenly 
transferred from the condensate |2), initially in its ground 
state in the trap, to the condensate |1). Due to the differ- 
ent magnetic moment, the condensate |1) experiences a 
trapping potential having the axial and radial frequencies 
lowered by a factor \[2 and the minimum shifted radially 
by a constant value. Thus, the trapping potential V* are 
different and can be written as 



Vt(x, y, z) = -lj 2 ±1 [x 2 + (y + y Q ) 2 + \ 2 z 2 ] (2) 
F 2 ( W ) = ^i 2 [z 2 +y 2 + AV] . (3) 



From now on we use the same parameters as in the ex- 
periments of Ref. H, that is: scattering lengths 022 = 
Oi2 — 98ao and an ~ 94.8ao, trapping frequencies 
lo ±2 = 2tt x 164.5 Hz, oj ±1 = u} ±2 /V2 = 2n x 116.3 
Hz, and displacement yo — g/u)j_ 2 — 9.2/xm. The 
two condensates have the same asymmetry parameter 
A = oj z /oj± ~ 0.0766. For the total number of atoms 
we use N = Ni + N 2 = 1.5 X 10 5 with N ± = 0.137V. 

In order to numerically solve the GP equations ([!]) one 
should map the two order parameters on a fully three- 
dimensional grid of points. This should be large enough 
to describe the oscillations of each condensate in the trap 
(center-of-mass motion plus internal excitations) as well 
as their expansion when the trap is switched off. This ap- 
proach is computationally heavy. We prefer, at this stage, 
to solve the GP equations ([!]) within different approxi- 
mate schemes, depending on the type of process we want 
to describe. We distinguish the following three cases. 

Case A: dynamics of a condensate in the trap ignoring 
its interaction with the other. In this case, the center-of- 
mass (CM) motion is trivial: the CM of the condensate 
\2) remains at rest, while the one of |1) oscillates freely 
around yo at the frequency Wxi- The internal motion of 
each condensate originates from the fact that its start- 
ing configuration, after transferring atoms from |2) to 
|1), is not the ground state configuration. In particu- 
lar, the state |2) feels the same trap but has less atoms: 
its starting shape is wider than the equilibrium one. Its 
internal motion can be predicted by exactly solving the 
corresponding GP equation, without coupling with |1), 



i.e., with cii2 = 0. In this case, the equation has axial 
symmetry and we can use the numerical code developed 
by two of us in Ref. @. 

Case B: effects of the mutual interaction in the radial 
motion. In order to point out the basic effects of the 
interaction occurring between the two condensates when 
they collide, one can decouple the slow axial motion from 
the GP equation. In this way, the simulation is reduced 
to a two-dimensional problem in the xy-plane, which is 
easier to perform. We solve the coupled GP equations 
by using the Split-Step Fourier method pT| , with a Fast 
Fourier Transform algorithm ]l2|] , by using cylindrical 
coordinates and dimensionless variables, as explained in 
the Appendix. We will justify the decoupling of the axial 
motion later, when discussing the low energy collective 
modes of these condensate. 

Case C: dynamics of the free expansion. In the ex- 
periments, the CM positions and the aspect ratios of the 
condensates are measured after few milliseconds of ex- 
pansion, with the trap switched off. We simulate the 
expansion, starting from the configuration obtained in 
case A or B at a given time. The CM of the two con- 
densates simply undergo a free fall under gravity. The 
shape of each condensate changes during the expansion. 
We simulate the expansion of the condensate |2) by using 
a scaling law valid in the Thomas-Fermi (large N) limit 
We instead ignore the mutual interaction of the 
condensates in this process; even if they overlap after 
some time (see Section [y]), their density is much lower 
than the initial one, due to the fast radial expansion, and 
hence the mean-field interaction is also much weaker. 



III. CENTER-OF-MASS MOTION IN THE TRAP 

We start from an initial configuration corresponding to 
the lowest stationary solution of Eq. (|I|) , with N atoms 
in the |2) condensate. This can be obtained by iterating 
in imaginary time, starting from a trial wave function, 
as done in [ p^[ . Then, at an arbitrary time t — 0, one 
builds two condensates |1) and |2), with the same density 
profile of the starting one, but with different number of 
atoms, Ni and N 2 . 

If they do not interact (ai 2 = 0), the CM of condensate 
|2) remains at rest, while the one of the condensate |1) 
starts oscillating around its equilibrium position y = — yo 
with frequency oj±i, as shown in Fig. [l]a. 

In order to include the effects of the mutual interac- 
tion, we make a two-dimensional simulation as explained 
in case B of the previous Section, by ignoring the ax- 
ial motion and solving the two coupled GP equations 
in the xy-plane. Thus, the initial condensate is an infi- 
nite cylinder along z while the normalization of the two- 
dimensional order parameter in the sy-plane is chosen 
in such a way to preserve the profile in the radial direc- 
tion, when passing from 3 to 2 dimensions, as explained 
in the Appendix. The results are shown in Fig. [IJ). One 
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sees that now both condensates move. The motion of 
the condensate |2) originates from the mutual repulsion 
with the condensate |1), which acts both at the begin- 
ning, when |1) is formed and starts moving away from 
1 2), and when they periodically collide. The motion of 
|1) is still dominated by large oscillations around —yo, 
but their frequency u>i turns out to be larger than the 
non-interacting frequency u)±\ by about 5.4%. Further- 
more, the oscillations appear to be damped. We notice 
that, though the minimum separation between the two 
center-of-mass positions slowly grows in time, it remains 
always smaller than the sum of the two condensate radii, 
at least in time interval here considered, thus ensuring 
that the two condensates indeed collide. 

The shift in frequency of the CM oscillations is one 
of the clear signatures of the the mutual interaction. It 
originates from the fact that the condensate |1) moves in 
an effective trap which is anharmonic. In fact, the pres- 
ence of the condensate |2) at one side of the harmonic 
potential in which |1) moves, provide an extra repulsion 
with a consequent up-shift of frequency. This shift only 
occurs if the two condensates overlap periodically. Since 
it comes from the cross mean-field term in the GP equa- 
tion, it depends significantly on A^. 

Indeed, as a first approach we can consider the ID 
problem of one single classical particle moving in the ef- 
fective potential V e (y) given by the harmonic oscillator 
plus the repulsion of the |2) atoms: 

K(2/) = ^i l2 ; 2 -|- 521 n max|o,l-^^| (4) 

where <?2i = ^fi 2 a2i/m, while R±2 and no are the ra- 
dius and the central density of the 1 2) condensate in the 
Thomas- Fermi (large N) limit, respectively |Q. To ob- 
tain the oscillation period, we carry out two integrals 
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where yi indicates the classical inversion point, equal to 
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The period depends on the number of atoms only through 
the radius R±2 that scales as N^ 5 . In Fig. || we see that 
the frequency shift predicted by Eq. (^|) is nicely close to 
the GP results. 



One can better visualize the effects of the interactions 
in the collisions by plotting the various energy terms con- 
tributing to the total energy of the system, as in Fig. ||. 
We checked that the total energy is accurately conserved 
during the simulation. The energies of each condensate 
(kinetic, mean-field, trapping) vary with time in a char- 
acteristic way: they change only when the condensates 
collide. This happens about every 8 ms. During the short 
collision time, when they partially overlap, the mutual in- 
teraction energy, i.e., the term containing the scattering 
length ax2, is significantly non-zero. In the upper part of 
the same figure we redraw the CM position of the con- 
densate 1 2), with a vertical scale larger than in Fig. |l|. On 
notices that the slope of these oscillations changes when- 
ever the two condensates collide, the precise change (sign 
and amplitude) depending on the phase of the oscillation 
when each collision occurs. 

The behavior of the energy vs. time is strictly related 
to the damping of the CM oscillations of condensate |1), 
shown in Fig. |l|b. First, we remind that in the GP theory 
the damping is not due to dissipative processes, since 
the total energy is conserved. What happens, instead, 
is that the kinetic energy initially associated with the 
CM motion is eventually shared among many degrees of 
freedom. This redistribution of energy is caused by the 
periodic collisions. One of these degrees of freedom is 
the CM motion of the condensate |2). But the internal 
excited states of each condensate are also expected to 
play an important role, as discussed in the next section. 



IV. CONDENSATE DEFORMATIONS AND 
ASPECT RATIOS IN THE TRAP 

Let us consider the internal motion of the condensate 
1 2). Just after the transfer of N\ atoms into |1), its den- 
sity distribution has the same profile as the ground state 
of N atoms, but with 13% of atoms missing. Thus, it 
starts oscillating around the new equilibrium configura- 
tion, the confining potential being unchanged. This sit- 
uation is very similar to that studied in Ref. Jl6[ , where 
similar oscillations were produced by a sudden change of 
scattering length. 

If one ignores the interaction with the condensate |1), 
the motion of the condensate |2) can be predicted by 
solving the GP equation ([!]) with i = 2 and 012 = 0. This 
corresponds to case A in Section 0. In Fig. |, we show 
the resulting aspect ratio, [(y 2 )/(z "j 1 / 2 (solid line). The 
curve is clearly a superpositions of two, fast and slow, 
modes. 

To understand which excited states partecipate in the 
curve in Fig. ^, we notices that, since the new equilib- 
rium state is not far from the starting configuration, the 
amplitude of the induced oscillations is relatively small. 
Therefore, only the lowest energy modes are excited and, 
among them, the ones with no angular momentum. For 
an elliptic trap, where only the z-component of the an- 
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gular momentum is conserved, two m = states are ex- 
pected to be excited. The one with lowest frequency 
corresponds to in-phase oscillations of the width along 
x and y and out-of-phase along z. The one with highest 
frequency is an in-phase compressional mode along all di- 
rections (breathing mode). In the Thomas- Fermi (TF), 
large N, regime and in the linear limit (small amplitude), 
their frequency is jl7j 

uj t = V<7t w -L2 (7) 

with 

q T = 2 + (3/2)A 2 =f (1/2) ^9A 4 - 16A 2 + 16 . (8) 

When the condensate is also strongly elongated (Acl) 
the two frequencies become 

uj + ~ 2 uj±_ 2 ; w_ ~ A cj ±2 • (9) 

In this limit the two frequencies are quite different, and 
the axial and radial collective excitations are almost de- 
coupled. The radial width is characterized by small os- 
cillations of frequency w_ superimposed to a wider os- 
cillation of frequency u>+, and vice- versa for the axial 
width pi . Using the parameters for the condensate |2) 
in Eqs (@)-(||), one finds values which differ less than 
0.04% from their asymptotic A — ► values in Eq. (^|). 

The TF approximation corresponds to neglecting the 
laplacian of the density (i.e., the quantum pressure) in 
the GP equation. In this case, the density profile of 
the condensate in each direction is an inverted parabola 
which vanishes at the classical turning point. In the ra- 
dial and axial directions, these points, Rj_ and Z, are 
given by the condition /i = mui± 2 R'j_/2 = mw 2 2 Z 2 /2, 
where /i is the chemical potential fixed by the normaliza- 
tion of the order parameter to A^- A nice feature of the 
TF regime is that both the free expansion and the lowest 
collective modes exhibit an exact scaling behavior: the 
parabolic shape is preserved and the widths R± and Z 
scale as R ± (t) = R±(Q)b±(t) and Z(t) = Z(0)b z (t), The 
GP equation for the condensate |2) then transforms into 
two coupled equations for the scaling parameters b±(t) 
and b z (t): 

b ± + Jib ± -*.=o (10) 

b z +u;%-^ = 0. (11) 

The numerical solution of these differential equations is 
straitghforward. The result for the evolution of the as- 
pect ratio of condensate |2), R±/Z, is shown as a dashed 
line in Fig. || and compared with the one obtained with 
the exact solution of the GP equation (solid line) for the 
same condensate, [(y 2 ) / (z 2 )] 1 / 2 . The shape of the two 
curves is very similar as a consequence of the fact that 



the parameter Na/a^o, with ay lo — [h/ '(mujho] 1 ^ 2 , is in- 
deed much greater than 1, as required to ensure the va- 
lidity of the TF approximation. They differ only for a 
small vertical shift due to the different definitions of the 
aspect ratio; this shift would vanish only in the N — > oo 
limit H). 

The aspect ratio in Fig. |^ ignores the mutual inter- 
action between the two condensates in the periodic col- 
lisions. Since this interaction breaks the axial symme- 
try in the motion of |2), one should solve the coupled 
GP equation in full tree-dimensions. Here we proceed 
differently. We observe that the axial motion is much 
slower than both the time interval between two collisions 
and the duration of a single collision. We thus assume 
that the axial motion does not respond significantly to 
the mutual interaction in the time scale here considered. 
Vice-versa, the collision are expected to significantly af- 
fect the radial motion. We treat the radial motion by 
solving numerically the GP equation in the xy-plane, in- 
cluding the mutual interaction (case B of Section |TJ) . In 
this case, the collisions with the |1) condensate make the 
radial motion of |2) rather complex. Assuming the ax- 
ial width to be the same as without collisions, and using 
the radial width obtained by the 2D integration of the 
coupled GP equations, one can calculate the new aspect 
ratio. This is given in Fig. |^. One notices that, similarly 
to what we have found for the center-of-mass motion in 
Fig. |l|b, the shape of the oscillation change at each col- 
lision, every 8 ms. The evolution of the aspect ratio has 
a complex structure and it is not a simple superposition 
of the two oscillations w_ and lu + of the non-interacting 
case. The overall amplitude is enhanced. This reflects 
the fact that the energy of the CM oscillation of |1) is 
transferred, not only to the CM motion of |2), but also 
to the internal degrees of freedom. This again, may ex- 
plainpart of the damping of the CM oscillations of |1) in 
Fig. 0b. 

Looking at Fig. [| one sees also that the interaction be- 
tween the two condensates produces an initial delay for 
the onset of collective exitations, compared to the non- 
interacting case. This originates from the fact that the 
two condensates are initially created with the same den- 
sity profile, and the scattering lenghts CI22 and an have 
the same value. Therefore the |2) condensate "realizes" 
to be out-of-equilibrium only when the other condensate 
has moved away, i.e. after a time of the order of 1 -7- 2 
ms. The same effect can be seen in Fig. || for the center- 
of-mass motion. 



V. FREE EXPANSION AND COMPARISON 
WITH EXPERIMENTS 

In the experiments of Ref. Q the center-of-mass mo- 
tion and the aspect ratios are measured after releasing 
the trap and letting the condensate expand under the 
effect of gravity only. 



4 



In order to compare our theoretical predictions with 
the experimental data we have to include the effects of the 
expansion. For the center-of-mass motion this is trivial. 
One has only to take care of the fact that the switch-off 
of the trapping potential is not istantaneous and can pro- 
duce a magnetic gradient which affects the initial veloc- 
ity of the two condensates. The two velocities have been 
measured experimentally, giving v\ v = 0.7±0.1 cm/s and 
V2 y = 1.4 ±0.1 cm/s. Due to this difference, the two con- 
densates fall at different speed and they eventually over- 
lap for some time during the expansion. As already said 
in Section || (case C), we ignore the mutual mean- field 
interaction during this process. If the overlap does not 
occur too early, the condensates have enough time to ex- 
pand fast in the radial direction, the one of the stronger 
initial confinement, and their density drops quickly, so 
that the cross mean-field energy becomes negligible. 

Our results for the CM position of the two condensates 
are shown in Fig. ^. In the upper part (a), we compare 
theory and experiment in a simple case, in which all the 
atoms are initially transferred into the condensate |1), 
which then oscillates alone in its new trap. The oscilla- 
tion is undamped and its frequency is simply the trapping 
frequency u>±i — 2n x 116.3 Hz. This picture, in which 
the mean-field interactions does not enter at all, has been 
used to fix the small delay caused by the non istantaneous 
switch-off of the trapping potential. This has been done 
by introducing the delay as a fitting parameter in the 
theoretical curve. Its value is ~ 0.6 ms. 

In the lower part of Figs. |^, we compare theory and 
experiment for the case of the two colliding condensates. 
The CM dynamics of the two condensates in the trap has 
been discussed in Section III. Starting from each point 
in Fig. [[Jd, we have then added the effect of the fall in 
the gravity field. At the time of the observation (after 
29.3 ms), the condensate |1) is almost always above the 
condensate |2), which means that they crossed during the 
expansion. 

There is an overall qualitative agreement between the- 
ory and experiments. The CM of the condensate |2) 
moves, due to the interaction with the other conden- 
sate. The amplitude of the oscillation is in agreement 
with theory, even though more precise measurements are 
needed to make the comparison more quantitative. The 
CM of the condensate |1) exhibits a damped oscillation. 
Its frequency is up-shifted compared to the case with 
no collisions (upper part of the same figure). The mea- 
sured frequency in the case of colliding condensates is 
123.9(3) Hz, with a shift about 6.4(3)%, which is close to 
the theoretical shift, 5.4%. Furthermore, the oscillations 
appear now damped with an exponential decay time of 
about 60 ms. Our theoretical prediction for this decay 
time is about 170 ms. As already said in the previous 
sections, the predicted damping originates from inelastic 
processes in which the energy of the CM motion is trans- 
ferred to other degrees of freedom. The fact that our 
theory underestimate the decay time is reasonable if one 
consider that the theoretical curve in Fig. H only includes 



the radial degrees of freedom (we solve the coupled GP 
equations in the xy-plane). The remaining contribution 
to the damping can naturally arise from the excitations 
of axial modes, which are ignored in the present analysis, 
namely, the possible distortions of the shape of the two 
condensates along z, during each collision. A source of 
damping can also be the occurrence of elastic scattering 
processes between atoms of the two overlapping conden- 
sates, as the one discussed in ppf . These processes are 
not included in the GP mean-field theory. 

In order to compare theory and experiments for the 
measured aspect ratio, one has to include the effect of 
the mean-field interaction in the expansion. For the con- 
densate 1 2) one can safely use the TF approximation. In 
practice, one has to solve the coupled equations ( |il"|) for 
the scaling parameters b± and b z , after dropping the term 
linear in &, which represents the trapping potential. The 
initial conditions for b(0) and 6(0) are fixed by the con- 
figurations of the condensate |2) in the trap at a given 
time t, as calculated in Section IV (see Fig. ||). In Fig. 
^ we plot our results, as a function of t, compared with 
the experimental data. The oscillations are rather com- 
plex, as a result of the periodic collisions between the 
condensate. The average amplitude is in agreement with 
the observed data. The slow axial modulation is visible 
and the faster radial motion has sudden changes of am- 
plitude. Considering that in this comparison there are 
no fitting parameters, and the theory contains some im- 
portant approximations, the overall agreement is good, 
and provides a further evidence of the mutual interaction 
between the two colliding condensates. 

Finally, let us comment on the aspect ratio of the con- 
densate |1). This quantity is expected to follow an even 
more complex evolution, due to the fact that the initial 
configuration is very far from equilibrium: the trap fre- 
quencies are changed by a \/2 factor, the condensate is 
initially displaced from the minimum of its trap and has 
a shape originally build for N atoms instead of N±. In 
the TF limit these reduction of the trapping frequencies 
and of the number of atoms would be equivalent to an 
"effective" sudden change of the trapping frequencies, of 
about 96%, if N were kept fixed. Therefore, due to the 
non-linearity of the GP equation, the resulting evolution 
might be chaotic, as predicted in Ref. Even though 
the TF approximation is rather poor for this conden- 
sate, which has about 20000 atoms, we solved the TF 
equations for the evolution of R±/Z finding rather irreg- 
ular oscillations between the extreme values 1 -j- 15. The 
comparison with experimental data is thus difficult, even 
because aspect ratios greater than about 5 are not easy 
to be measured with the finite resolution of the imag- 
ing devices of Ref. [Q . In view of the characterization of 
the nonlinear (possibly chaotic) dynamics of these con- 
densates, we think that a deeper quantitative analysis of 
this problem is worth doing in the next future by solving 
the GP equation in full 3D. 
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APPENDIX A: SOLVING THE 
TIME-DEPENDENT GP EQUATION 

The GP equations in Eq. (Q) can be conveniently 
rewritten in terms of dimensionless quantities, by rescal- 
ing all time and space variables by lo±_ = lo±_2 and 
a± = ^Jh/muj\_ respectively. In fact, by defining 

^ = x/a±; £, 2 =y/a±; £ 3 = z/a ± (Al) 



8t 



.d(f)2 



3 = 1,2 



-V 2 + £ 2 + (6 - 6o) 2 + A 2 £ U 2 j\<f>jf 

3 = 1,2 



(A8) 

and we solve them by using the Split-Step Fourier method 
iPH , with a Fast Fourier Transform algorithm |l2j, and 
with the normalization condition (i = 1,2) 



d 2 Z \MZ)\ 2 = 1 



(A9) 



w x t/2; 




8irNj ciij /a± 



we can write Eq. (Q) as follows 



1 dr 



3 = 1,2 



(A2) 



ULr) (A3) 



where the magnetic trapping potential are given by (A 



^(0 = 2 [ei 2 + (6 + 6o) 2 + A 2 C 2 ] 

v^) = £ 2 + £ 2 2 + a 2 £: 



3) 



(A4) 
(A5) 



and the wave functions ipi satisfy the normalization con- 
dition (i = 1,2) 



rf 3 £ \M0\ 2 = i. 



(A6) 



One of the scheme we use (see Section [Oj) consists in 
reducing the dynamics to the sy-plane, by neglecting the 
axial motion. This corresponds to solve the dynamics of 
the trapped oscillations in the limit A = 0, i.e., assuming 
a uniform wave function along the axial direction. We 
write ip = -<4<^(£i, £2)1 and we fix A in order to preserve 
the width of the condensate profile in the radial direction, 
when passing from 3 to 2 dimensions (in the TF limit): 



A _ (15A) 2 / 5 ^Na 22 ^ 
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(A7) 



Then the GP equations (|A3|) become 
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FIG. 1. Center-of-mass oscillations of the |1) (solid line) 
and |2) (dashed line) condensates in the trap: (a) ignoring the 
mutual interaction, (b) including the mutual interaction. The 
dotted line corresponds to the equilibrium position y — — yo 
of the |1) condensate. Time and lengths are given in units of 
and a±2 = [h/(muj±2)] 1 ^ 2 , respectively. 
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FIG. 2. Frequency shift of the center-of-mass oscillation of 
the |1) condensate as a function of the number of atoms N2 
remaining in |2). GP results (points) are compared with the 
classical prediction of Eq. (^) (solid line) . 




FIG. 3. Plot of the center-of-mass oscillation of |2) con- 
densate (top) and of the various energy terms contributing to 
the total energy of the system (bottom), vs. time. The total 
energy (dotted line) is divided as follows: energy of the |2) 
condensate (kinetic + potential + mean-field, dashed), energy 
of the |1) condensate (dot-dashed), and interaction energy be- 
tween the two condensates (solid line). 
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FIG. 4. Evolution of aspect ratio for the |2) condensate in 
the trap, as a function of time (in units of w7_ 2 ) when the 
interaction with the |1) condensate is neglected. The solid 
line is the GP prediction for [{j/ 2 )/(« 2 )] 1/ ^ 2 , while the dashed 
line is R±/Z in the Thomas- Fermi limit. 




FIG. 5. Evolution of aspect ratio for the |2) condensate in 
the trap, as a function of time (in units of wjj). The P er i°dic 
interaction with the |1) condensate is included in the case of 
the solid line and neglected for the dashed line. 
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FIG. 6. Center-of-mass oscillations as a function of the 
trapped evolution time (in units of wT^), after t exp = 29.3 ms 
of ballistic expansion: a) CM oscillations of the |1) conden- 
sate alone; b) CM oscillations of both condensates, subject 
to mutual interaction. Theory (solid lines) is compared with 
experiments (points). 
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FIG. 7. Aspect ratio for the |2) condensate, 



func- 



tion of the trapped evolution time (in units of w^ 2 )' after 
texp = 29.3 ms of ballistic expansion. The experimental points 
are compared with the theoretical predictions for the case of 
two interacting condensates. 
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